Instability of Chaplygin gas trajectories in unified dark matter models 
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In the past few years, the Chaplygin gas (CG) has been considered as an appealing candidate 
for unifying dark matter and dark energy into a single substance. This picture had to face several 
problems when trying to fit its predictions with cosmological observations. We point out another 
potential problem of this model, showing that, when the CG is described through a classical scalar 
field, the corresponding trajectories are strongly unstable. This implies that an extreme fine tun- 
ing must apply to the initial values of the field in order to end up today with values allowed by 
observational data. 



Since the observations of distant type la Super- 
novae pointed towards an accelerated expansion of 
the Universe 0, 1 a great effort has been done to 
provide a reasonable, physically motivated respon- 
sible for this speeding up. Apparently, a dark en- 
ergy component adds up to the dark matter, coming 
to dominate the total energy budget at very recent 
times. The remarkable new feature of dark energy 
is that it appears to violate the strong energy con- 
dition 3]; understanding the nature of dark energy 
is probably one of the most important open prob- 
lems in modern physics. Besides the possibility of 
a quintessence scalar field 4] , modifications of grav- 
ity 0, or an uncanceled cosmological constant 0, 
an interesting alternative class of dark energy mod- 
els is that of the so-called Chaplygin gas 7]- 23]. 
In its simplest formulation, the Chaplygin gas is a 
perfect fluid with equation of state p = —A/p, A 
being a positive constant with the dimensions of an 
energy squared. A generalized form of the Chaply- 
gin gas, containing an additional free parameter, has 
also been studied in detail (see, e.g., 10, 14]). 
While first introduced in a hydrodynamics context 
[2^ . the Chaplygin gas has recently raised growing 
interest in particle physics, thanks to its connection 
with string theory [1 HI HHHIli HI |^ and 
because it is the only fiuid known to admit a super- 
symmetric generalization jsj]- 

From a cosmological point of view, the striking fea- 
ture of the Chaplygin Gas (hereafter, CG) is that it 
allows for a very elegant unification of dark energy 
(DE) and dark matter (DM), since its equation of 
state interpolates between a dust dominated phase 
at early times and a de Sitter phase at late times. 
Despite the appeal of these intriguing features, sev- 
eral fatal drawbacks turned out to affect CG mod- 



els. There are, actually, two different approaches 
to the cosmological role of the CG: one can con- 
sider it as a unified dark matter component, so that 
^DM = and the CG plays both the roles of dark 
matter and of a cosmological constant, at different 
epochs (see, e.g., Ref. !32|); or one can look at the 
CG as a Dark Energy candidate, which adds to the 
standard DM component, as well as to baryons and 
radiation [l^ Il9|. In both cases, the CG is unable 
to match several observational data. 
In Ref. [i^l it was pointed out that perturbations 
in the CG should affect the formation of structures, 
producing oscillations or exponential blowups in the 
matter power-spectrum which are inconsistent with 
observations; the analysis ruled out the CG as a uni- 
fied dark matter candidate {^Idm = 0) at 99.99% 
level. 

The impact on cosmology has also been investigated 
in the context of Supernovae data 0, 33, 34 
I35I liif . showing initially a certain degree of con- 
sistency with data for a class of CG models; Cos- 
mic Microwave Background (CMB) measurements, 
though, provided stronger constraints to CG cos- 
mologies [H El EH. While Ref. jl^ found that 
the joint analysis of type la Supernovae and CMB 
data only allows for a CG gas which is today indis- 
tinguishable from a cosmological constant, the anal- 
ysis of Ref. .18] rules out the CG as a Dark Energy 
candidate at 99.99% confidence level, and also indi- 
cates that the CG as a unified model of dark matter 
and dark energy is strongly disfavoured by the latest 
CMB data. 

In this paper, we point out another possible prob- 
lem related to the dynamical behavior of the CG 
gas. We focus on the unified dark matter cosmo- 
logical model and we analyze the trajectories of the 
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classical scalar field which is meant to correspond 
to the CG hydrodynamical representation, search- 
ing for attractor solutions. For a flat universe, the 
parameters of the CG can be easily related to the 
expansion rate today, Hq , and to the present cosmic 
equation of state, wq. The latter quantity is con- 
strained by the recent estimates of the present de- 
celeration parameter qq since wo = (2^0 ^ l)/3 
; we thus require wq —0.7 today. We find that, in 
order to finish up at the present epoch with such an 
equation of state wq (different than —1) and a plau- 
sible value of the expansion rate, a strong fine-tuning 
must be applied on the initial conditions of the scalar 
field. In other words, the trajectories which repro- 
duce the observed parameters are strongly unstable 
under variations of the initial conditions on the field. 
The only attractor of the dynamical system under 
study is indeed the de Sitter Universe, to which the 
model converges at late times. Our conclusions are 
complementary to the perturbation analysis of |3^ . 
showing that, already at the background level, the 
unified dark matter model shows serious instabili- 
ties. 

We will undertake our analysis in a fiat 
Friedmann-Robertson- Walker (FRW) cosmology; 
since we are mostly interested in the dynamics at 
recent times, radiation is not included in our treat- 
ment. From the energy conservation and the CG 
equation of state, it follows that 
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a being the scale factor and B an integration con- 
stant; the ratio B/A characterizes the transition be- 
tween the matter-like behavior and the cosmological- 
constant. 

In the unified dark matter model, the expansion rate 
and the cosmic equation of state today are, respec- 
tively, = \/A + B and wq = -(1 + B/A)-^ . 
We see that the constants A and B can be deter- 
mined once Hq and wq are known. Our aim is to 
provide a description of the Chaplygin gas through 
a minimally- coupled, classical scalar field. In |7|.p^ 
this description is simply built from the analogy of 
the energy density and pressure of the CG in Eq. 
with the analogous quantitites for a scalar field (j) 
evolving in a potential V{4>), i.e. 



= 0V2 + 
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where a dot denotes differentiation w.r.t. the cosmic 
time. 



Namely, by requiring that = pchap and = 
PC hap, we can relate the cosmic scale factor to the ki- 
netic and potential energy a scalar field should have 
in order to reproduce the Chaplygin fiuid dynamics: 
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By making use of the Friedmann equation and of 
Eq. one can infer d(/)/da, which can be solved to 
give 0(0); The authors of Ref. 3 gi'^s 
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There is a crucial assumption in Eq. (@J: the inte- 
gration constant has been choosen so that the value 
of the scalar field at some (arbitrary) time ti is fixed; 
in particular, setting the initial conditions at a^, one 
should have 
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The initial "velocity" (^i is fixed too, by Eq. (jSJ. 
Substituting a(0) in Via), one has 



y(0) = \^fA f cosh30 + 

2 \ cosh3(p 



(6) 



As already noticed in Ref. [23] , a scalar field evolving 
in such a potential will produce a cosmological evo- 
lution coinciding with that of the Chaplygin gas only 
if the initial conditions are appropriately handled to 
satisfy the relation (/)| = 4(l/^((/)i) — A); indeed, the 
potential JHIl has been built under this assumption. 
The evolution of a minimally-coupled scalar field fol- 
lows the Klein-Gordon equation of motion: 



-3H<I>-V4 



(7) 



As we said, in order for the solution of this equa- 
tion to reproduce the energy density and pressure 
Q), one has to exactly fix the initial values of (f>,(f>- 
This is quite intuitive, and it can also be easily un- 
derstood by looking at the shape of the potential 
lO, which is extremely flat for small values of the 
field, increasing exponentially elsewhere; as an ex- 
ample, a scalar field initially set on the minimum 
of this potential with vanishing velocity, would ex- 
hibit no dynamics at all, and it would be indistin- 
guishable from a cosmological constant at any cos- 
mic epoch. Of course, this is not what we are ex- 
pecting from a plausible scalar-field representation 
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of the Chaplygin fluid; first, a unified model of dark 
matter needs at least an epoch when the equation of 
state Wij, = p^/prj, is zero, so that the energy density 
drops like a matter component. This phase should 
end up in a de Sitter phase, when the energy den- 
sity of the scalar field is constant (as for a >> 1 in 
(P) . Finally, since the cosmic equation of state today 
is kinematically related to the deceleration parame- 
ter q = d/oH'^, the most recent measurements give 
Wo ~ 0.7 I33. It follows that the intermediate stage 
of the CG evolution should be characterized by a 
decreasing equation of state, shifting smoothly from 
to —1. The point we focus on here is how much 
fine-tuning, if any, is required to do that. In other 
words, is there any freedom in setting the initial val- 
ues of the field if one has to reproduce the scalings 
in Eq. (Q? The question may be turned into the 
phase-space language, by asking whether the Chap- 
lygin fluid trajectory is an attractor for the system 
O. As we said, we can always provide a scalar- 
field description of the Chaplygin gas, by appropri- 
ately setting its initial conditions and building up 
the corresponding potential. Such a field will obey 
the Klein-Gordon equation, and, in addition, p^,Pcj, 
will scale as in Eq. at any time. The scalar field 
model obtained with this particular choice of initial 
conditions would perfectly reproduce the dynamics 
of the Chaplygin gas. We will refer to (jichapit) as 
the solution of Eq. |(7|) with potential ((HJ, which sat- 
isfies the required initial conditions. 

As for the late time behavior of Q , it is easily seen 
that, whatever initial conditions we choose, the sys- 
tem will end up in its global attractor, the de Sitter 
stable node (see [s^) with (j),<j) = 0, H = \fA (note 
that the de Sitter node is the only critical point of 
the dynamical system under investigation). In par- 
ticular, this is true for the solution (/)chap(t). How- 
ever, we are interested in finding whether the tra- 
jectory (jichap itself is an attractor for the problem 
on hand, in the sense that many possible solutions of 
the system converge to a common evolutionary track 
represented by the 4>chap trajectory. We find that 
this is not the case; while 4'Chap requires fine-tuned 
initial values to reproduce the Chaplygin gas cosmol- 
ogy, slightly different initial conditions on 4>Chap may 
result in completely different trajectories, which may 
even be unable to reproduce the matter behavior at 
all. The general solutions of this dynamical system 
turn out to behave much differently than the "track- 
ing solutions" of • Following the notation of i 
we may anticipate our result, saying that the system 
does not admit tracking solutions because one of the 
two necessary conditions on F = VV^^/V^ is vio- 



lated, namely F = 1 for almost any plausible initial 
field value. One should be careful, however, in using 
the formalism of js^ in this case, in that the Chap- 
lygin gas is supposed to be the only component of 
the cosmic fluid, and there is no other background 
component which may drive its trajectories (having 
neglected the role of radiation in our analysis). For 
this reason we found it useful to directly check for 
the stability of the 4>Chap solutions with a numerical 
integration of the system. 

Our plots refer to four different choices of initial 
conditions, set at redshift z ^ 10^. The solid lines 
refer to the 4>chap trajectories. The present values of 
Hq and wq unambiguously determine A and B for a 
flat universe and the present values of 4'Chap, 4'Chap- 
We set h = 0.7 and wq = -0.7. Together with the 
4chap trajectories. Figs. Q and El show, respectively, 
the energy density of the scalar fleld and the cor- 
responding equation of state, when the fleld starts 
with the following initial conditions: 

a) 4>i = 4>Chap,i X 1.1, (f>i^ 4>Chap,i (dotted lines); 

b) 4i = 4chap,i, 4i = 4chap,i X 10^ (dashed fines); 

c) 4i = 4chap,i X 0.9, 4'i. = 4chap.i X 10 (dot-dashcd 
lines) . 

Even a slight modification of the initial values with 
respect to 4'Chap can make the equation of state to- 
day completely different from the required value, re- 
sulting in trajectories which never join the 4chap 
curves, nor reproduce the CG behavior at any epoch. 
Furthermore, the conditions a), b), c) finish up at 
present with unacceptable values of the Hubble ex- 
pansion rate. Starting with a bigger amount of ki- 
netic energy, as in b) and c), an initial kination pre- 
cedes the phase in which the field behaves like mat- 
ter. In each of the plotted cases, there is a very sharp 
transition between the matter-like scaling {w = 0) 
and the A-like one {w = —1): because of this sharp- 
ness, values of the equation of state different than 
0,-1 turn out to be unnatural. This behavior is 
very different than in tracking quintessence solu- 
tions, where the only parameters which need to be 
adjusted are the amplitude of the potential and the 
current equation of state: in that case, once these 
present values are fixed, there is no dependence on 
the initial values of the quintessence field. 

Concluding, we have shown that, in order to 
have a unified dark matter model consistent with 
the most recent estimates of the cosmic equation 
of state, one is forced to introduce an extreme 
fine tuning on the initial values of the scalar field, 
because the problem does not admit tracking 
solutions; the de Sitter universe is the only global 
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FIG. 1: Energy densities for the (pchap solution and for 
different initial conditions, corresponding to the cases 
a),b),c) described in the text. 



FIG. 2: Equation of state for the (pchap solution and for 
different initial conditions, corresponding to the cases a), 
b), c) described in the text. 



attractor of the system. 

Although the potential © requires ad hoc initial 
conditions on 0, which allowed us to define (jichap, 
this is not a unique representation of the Chaplygin 
fluid. One may start with different values of 
and, by a similar approach, build a potential 
V{4)) (generally, different from O for each set of 
initial conditions, so to reproduce the behavior 
n In this way, there would be an infinite class of 
potentials, mapping the behavior of the Chaplygin 
fluid. However, our conclusions apply to any on 
this potentials, due to the peculiar sharpness of the 
matter-A transition of the CG fluid. 
Therefore, the CG as a unified dark matter candi- 
date, seems to be strongly disfavored from the point 
of view of its dynamics. 
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